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Final  report: 


Digital  Holographic  Interferometry  for  Airborne  Particle 

Characterization 


Matthew  J.  Berg 
Mississippi  State  University 
Department  of  Physics  &  Astronomy 


Abstract 

This  final  report  covers  the  research  activities  supported  by  grant  59999-EV-YIP.  In 
addition  to  launching  the  Pi’s  new  research  program  at  MSU,  three  graduate  and  five 
undergraduate  students  were  supported,  nine  peer-reviewed  journal  articles  were 
published,  and  20  presentations  at  professional  meetings  were  given.  Several  major 
research  discoveries  were  made  including  a  relationship  between  an  aerosol  particle’s 
hologram  and  its  extinction  cross  section,  and  a  computational  demonstration  that 
holographic  interferometry  can  resolve  aerosol  particle  size  evolution. 


Laboratory  infrastructure  and  educational  achievements 

The  first  major  accomplishment  of  this  work  was  to  establish  a  new  electromagnetic  laboratory 
for  the  PI  and  launch  his  research  program  in  holography  applied  to  aerosol  characterization. 
This  was  achieved  through  the  Young  Investigator  Program  (YIP)  and  supplements  from  the 
Defense  University  Research  Instrumentation  Program  (DURIP)  during  the  project’s  second 
term.  Equipment  obtained  with  this  support  includes,  e.g.,  UV  and  visible  diode-pumped 
Nd:YLF  Q-switch  lasers  from  Photonics  Industries,  an  optics  table,  optical  elements  with 
associated  hardware,  electronic  instrumentation,  large-format  CCD  sensor,  and  a  custom-made 
airtight  aerosol  cell,  see  Fig.  1.  As  the  laboratory  grew,  the  PI  recruited  graduate  students  to 
undertake  research  projects  related  to  the  YIP  proposal.  During  the  summer  term,  the  PI  also 
recruited  and  supported  undergraduate  students  to  work  with  the  graduate  students  through 
yearly  supplements  from  the  Undergraduate  Research  Apprenticeship  Program.  Several  of  these 
students  later  joined  the  graduate  program  at  MSU. 


Figure  1 :  Large-format  CCD  sensor  (a),  air-tight  aerosol  cell  (b),  and  the  new  lasers. 


Research  achievements 

The  research  objectives  of  the  original  YIP  proposal  were  pursued  through  multiple  laboratory 
and  computational  investigations  over  the  duration  of  the  grant.  Major  achievements  were 
reported  in  a  total  of  nine  peer-reviewed  journal  publications  and  described  in  presentations  at 
national  and  international  meetings  and  workshops,  including  twenty  presentations  in  all.  These 
presentations  were  highly  valuable  to  the  progression  of  the  work  as  input  and  advice  from  the 
Pi’s  colleagues  in  the  community  strongly  affected  how  various  technical  problems  were 
addressed.  In  what  follows,  a  brief  description  of  the  major  research  achievements  is  given.  More 
detailed  descriptions  can  be  found  in  the  published  journal  articles. 


Backscatter  digital  holography 

First  developed  was  a  method  for  digital  holographic  imaging  of  particles  using  their 
backscattered  light,  as  opposed  to  the  forward  scattered  light  pursued  previously.  This  effort  was 
motivated  by  the  potential  practical  applications  for  stand-off  aerosol  characterization.  In  short 
the  optical  arrangement,  which  was  a  Michelson-interferometer-like  setup,  used  a  beamsplitter  to 
direct  a  particle’s  back  scattered  light  onto  a  CCD  sensor.  Combination  with  the  retro-reflected 
incident  wave  across  the  sensor  produced  a  hologram,  from  which  a  particle  image  was  obtained 
following  the  same  computational  process  as  with  the  forward  scattering  technique,  see  Fig.  2. 
This  imaging  was,  overall,  successful  and  resulted  in  an  paper  published  in  Optics  Express. 
Another  motivation  for  this  work  was  to  determine  if  such  holograms  could  resolve  the  surface 
structure  of  opaque  particles.  Unfortunately,  this  question  could  not  be  answered  at  the  time 
since  the  resolution  of  the  images  was  too  coarse  to  make  such  an  assessment.  However,  the 
large-format  CCD  sensor  obtained  through  DURIP  support  following  these  experiments  could 
provide  the  needed  resolution  and  such  is  part  of  ongoing  work.  Also  notable  here  is  that  the 
work  was  partly  performed  by  a  former  URAP  student,  Mr.  Nicholas  Fowler,  who  enjoyed  his 
previous  URAP  experience  enough  to  return  for  a  second  summer  to  help  complete  the  work.  His 
contribution  was  acknowledged  in  the  publication. 
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Figure  2:  Diagram  (left)  of  the  experimental  arrangement  used  to  demonstrate  backscatter  digital 
holography  of  fixed  microparticles.  The  images  to  the  right  show  an  example  of  a  backscatter  contrast 
hologram  along  with  the  reconstructed  image.  The  particles  in  this  case  are  Aspergillus  flavus  spore  clusters. 


Holography  interferometr  y 

A  major  objective  of  this  project  was  to  investigate  the  possibility  of  holographic  interferometry 
to  infer  material-dependent  properties  of  aerosol  particles.  To  this  end,  a  computational  study 
was  first  undertaken  wherein  the  concept  was  tested  for  a  spherical  particle  undergoing 
expansion.  Initially,  it  was  believed  that  the  concept  would  readily  work  on  aerosol  particles  in 
the  laboratory  and  no  serious  computational  work  was  planned.  However,  following  multiple 


discussions  with  colleagues  at  professional  meetings,  the  feasibility  of  the  concept  became  less 
obvious,  and  so,  the  computational  study  was  justified. 

First,  to  describe  the  concept,  consider  the  illumination  of  a  particle  it  twice  and  allow  the  sensor 
to  record  a  double  exposure  consisting  of  the  two  contrast  holograms,  /1con  and  I2on.  For 
example,  these  illumination  stages  could  be  laser  pulses  timed  to  illuminate  a  particle  at  stages 
during  its  expansion  or  contraction.  To  better  understand  what  this  double-exposed  hologram 
represents,  suppose  that  the  particle  is  actually  in  exactly  the  same  physical  state  during  each  of 
these  illuminations,  i.e.,  has  not  changed  in  size.  In  this  case,  the  holograms  will  be  the  same  and 
application  of  the  Fresnel-Kirchhoff  (FK)  diffraction  integral  will  yield  two  overlapping, 
identical  images.  However,  if  the  particle’s  state  changes  between  the  pulses,  the  holograms  will 
then  differ.  Applying  the  FK  integral  in  this  case  will  yield  two  reconstructed  waves,  K1  and  K2 , 
that  differ  in  amplitude  and  phase.  Thus,  these  waves  will  interfere,  and  the  absolute  square 
\K1  +  K2\2  will  yield  a  particle  image  superimposed  with  interference  fringes.  The  fringes  are 
then  a  direct  consequence  of  the  change  in  the  particle’s  state  occurring  between  the  first  and 
second  illumination. 


To  test  this  concept,  an  expanding  spherical  particle  is  considered.  This  permits  the  use  of  Mie 
theory  to  generate  the  exact  scattered  fields  from  which  a  contrast  hologram  can  be  produced.  To 
do  this,  a  simulation  was  developed  in  the  C++  language  following  the  Mie  theory  derivation 
outlined  in  the  resulting  JQSRT  article.  Verification  of  the  code  is  done  by  comparing  the  far- 
field  scattered  intensity  patterns  for  a  variety  of  spheres  of  different  size  and  refractive  index  to 
that  produced  by  the  widely  used  BHMIE  code,  which  is  also  given  in  the  paper.  The  simulation 
considers  the  particle  with  two  sizes;  R1  and  R2  where  R2  >  R1  and  the  same  refractive  index 
m  —  1.33  +  0.05i.  Here,  the  scattered  electric  Esca  and  magnetic  Bsca  fields  are  computed  at 
each  pixel  in  the  hologram  sensor,  which  resides  in  the  particle’s  far-field  zone.  These  fields  are 
then  superposed  with  the  incident  plane  wave,  Einc  and  Binc.  From  this,  the  time-averaged 
Poynting  is  calculated  and  the  component  directed  into  the  sensor  is  assigned  to  each  pixel.  Thus, 
the  two  contrast  holograms  /1con  and  I2on  are  generated  corresponding  to  the  particle  sizes,  R1 
and  R2,  respectively.  These  holograms  are  then  superposed  to  yield  the  final  double-exposed 
hologram,  /tcc”tn .  To  this  the  FK  integral  is  applied  to  render  an  image  of  the  particle  as  \Ktot\  2. 
Figure  3  shows  how  a  change  in  the  particle  size,  i.e.,  A R,  appears  in  the  fringe  structure  in  the 
cross  section  taken  through  the  reconstructed. 


Figure  3:  Evolution  of  the  interference  fringe-structure  in  the  reconstructed  image  for  an  expanding  particle. 


Figure  4,  however,  shows  how  this  change  in  particle  size  is  expressed  in  the  two-dimensional, 
silhouette-like  reconstructed  images.  Here,  the  particle's  two  sizes  are  —  21 A  and  Rt  —  47 A, 
where  A  is  the  wavelength  of  the  incident  light.  The  first  two  columns  in  the  figure  display  the 
contrast  hologram  for  each  state,  /1con  and  I2on,  along  with  their  reconstructed  images,  shown 
below.  As  expected,  the  images  for  each  state  are  silhouettes  with  radii  Rt  and  R2.  For  the  state-2 
reconstruction,  the  dashed  white-outline  denotes  the  size  of  the  particle  in  state- 1  for  reference. 
Surrounding  each  of  these  silhouettes  are  faint  rings.  These  are  due  to  the  out-of-focus  twin 
image  that  is  always  present  with  the  in-line  geometry  used.  The  last  column  in  Fig.  4  shows  the 
double-exposed  hologram  7tc°tn  and  its  image  reconstruction,  |/ftol|2.  One  can  see  a  prominent, 
dark  inner-disk  surrounded  by  an  annular  region  with  lighter  fringes.  Here,  the  dashed  circles 
denote  the  Rt  and  R2  particle  sizes,  and  thus,  the  width  of  the  region  between  them  corresponds 
to  the  change  in  particle's  size  A R  =  R2  —  R1. 
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Figure  4:  Demonstration  of  holographic  interferometry  for  a  single  expanding  particle.  The  top  row  shows, 
respectively,  the  contrast  holograms  for  the  particle  in  state-1  I^on,  state-2  /£on,  and  the  double-exposed 
hologram  7tc0°tn.  The  bottom  row  shows  the  corresponding  reconstructed  particle-images  obtained  by  applying 
the  FK  integral.  The  particle's  refractive  index  is  m  —  1.33  +  0.05i,  and  in  state-1  has  size  R1  =  27  A  while 
state -2  has  size  R2  =  4771. 

As  Fig.  1  demonstrates,  the  change  in  size  of  the  expanding  particle  can  be  detennined  from 
either  the  number  of  interference  fringes  in  the  “overlap”  region  of  the  image  or  from  the  size  of 
this  region  itself.  Thus,  this  computational  effort  validates  the  basic  concept  proposed  for  this 
project  and  is  reported  in  more  detail  in  the  paper. 


Extinction  cross-section  extraction 

Extinction  is  a  ubiquitous  phenomenon  that  describes  the  attenuation  of  light  traversing  a 
medium  due  to  scattering  and  absorption,  and  the  ability  to  accurately  measure  the  extinction 
cross  section  Cext  is  important  in  many  applications,  including  aerosol  characterization.  An 
important  new  development  in  this  project  was  to  show  that  Cext  can  be  measured  by  recording 
the  interference  energy- flow  due  to  a  particle’s  scattered  light  interfering  with  the  incident  light 
across  a  detector’s  surface,  i.e.,  from  the  hologram. 


To  present  the  concept,  consider  a  single  spherical  particle  in  vacuum  illuminated  by  a  linearly 
polarized  plane  wave  traveling  in  vacuum  along  the  z-axis,  which  will  approximate  a  collimated 
laser  beam.  To  see  how  extinction  comes  about,  Poynting's  theorem  is  applied  to  the  volume 
defined  by  a  deformed  spherical  surface  consisting  of  component  surfaces  Sj  and  S2  shown  in 
Fig.  5  below. 


Figure  5:  Sketch  of  the  spherical  surface  5sph  and  composite  surface  51  U  S2  in  relation  to  the  particle  and 
incident  wave.  Here,  S2  is  coincident  with  the  surface  of  the  detector  used  to  record  the  particle’s  hologram. 


Using  the  divergence  theorem,  the  volume  integral  in  Poynting’s  theorem  can  be  transferred  into 
surface  integrals  over  and  S2,  giving 

l^abs  __jj  (S)-rd  a  +  J  (S)-zdaj,  Eq.  (1) 

where  I/Fabs  is  the  power  absorbed  by  the  particle  and  (S)  is  the  time-averaged  Poynting  vector. 
From  Eq.  (1),  the  following  quantities  are  defined: 

Idet  =  f  (Smc)-zda  and  Idet  —  f  (S)-zda,  Eq.  (2) 

Js2  s2 

which,  respectively,  represent  the  holograms  received  by  the  detector  without  and  with  the 
particle  present.  In  Eq.  (2),  Sinc  is  the  energy  flow  of  the  incident  wave.  Taking  the  difference  of 
these  two  measurements  with  Eq.  (1)  in  mind  gives 

7det  _  7det  =  ^abs  +  f  ^ext)  .  f  Eq. 

•'Si 

where  Sext  is  the  energy  flow  due  to  the  interference  between  the  incident  and  scattered  waves. 
Lastly,  we  define 

Cext  =  f(edet)  +  S(0det)  Eq.  (4) 


where 


and 


Eq.  (5) 
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Eq.  (6) 


where  /inc  is  the  incident  light  intensity  and  0det  is  the  angular  size  of  the  integration  surface  Sx 
on  the  detector. 


Figure  6:  Behavior  of  the  /  and  5  curves  of  Eqs.  (5)  and  (6)  as  a  function  of  ©det-  Plots  (a)-(c)  show  the 
curves  for  nonabsorbing  spheres  with  m  =  1.33  +  Ot  and  R=^.,  21.  and  51,  respectively,  whereas  the  sphere 
in  plot  (d)  is  absorbing  with  ?n  =  1.55  +  O.lt  and  R=51.  Also  shown  are  the  points  along  the  f  curve,  after 
which  its  error  in  approximating  Cext  drops  below  10%  (all  plots)  and  1%  in  plots  (c)  and  (d). 

What  remains  to  be  seen  is  exactly  how  this  extinction  behavior  is  related  to  holography.  To 
explain  this,  consider  a  Gabor-type,  or  in-line,  holographic  arrangement  intended  to  image  a 
single  particle.  Here,  a  laser  beam  directly  illuminates  a  position-sensitive  detector  facing  the 
oncoming  light,  i.e.,  what  is  shown  in  Fig.  5.  If  a  particle  is  introduced  into  the  beam,  the  light 
scattered  will  interfere  with  the  incident  light  and  produce  a  fringe  pattern  across  the  detector. 
This  interference  pattern  constitutes  a  digital  hologram  from  which  an  image  of  the  particle  can 
be  computationally  reconstructed  by  applying  the  FK  integral.  The  results  in  Fig.  6  show  that  by 
simply  integrating  the  contrast  hologram  from  the  forward-direction  out  one  can  get  an  estimate 
for  Cext.  The  accuracy  of  the  estimate  will  depend  on  the  particle  size  to  wavelength  ratio  and 
the  angular  size  of  the  detector,  but  as  shown  in  Fig.  6,  this  estimate  can  be  within  several- 
percent  error  from  the  true  cross  section  value.  Thus,  from  a  single  contrast-hologram 
measurement,  it  is  possible  to  extract  an  unambiguous  image  of  the  particle  simultaneous  with  a 


measurement  of  its  extinction  cross  section.  Further  detail  of  this  technique  are  given  in  paper 
published  in  Optics  Letters. 


Figure  7:  Spherical  void  electrodynamic  levitator  trap  (a)  and  high-voltage  power  supply  (b). 

To  test  this  concept  in  the  laboratory,  the  students  used  a  Spherical  Void  Electrodynamic 
Levitator  (SVEL)  to  collect  a  digital  hologram  of  a  single  glass  microbeads  suspended  in  the 
trap,  see  Fig.  7.  Then,  an  integral  of  the  hologram  from  the  point  corresponding  to  the  optical 
center  of  the  incident  beam  was  performed  to  yield  an  estimate  for  the  cross  section  of  the 
trapped  particle.  The  outcome  of  this  work  is  still  not  clear  due  to  the  difficulty  of  trapping  only 
a  single  particle  in  the  SVEL.  However,  work  with  this  experiment  is  ongoing  and  is  being 
extended  to  other  particle  types  that  may  be  easier  to  handle.  Once  the  experiment  can  be 
completed,  a  follow-up  to  the  Optics  Letter  published  on  the  concept  will  be  produced  that 
describes  the  experimental  tests. 


